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People's ability to think creatively is a primary means of techno- 
logical and cultural progress, yet the neural architecture of 
the highly creative brain remains largely undefined. Here, we 
employed a recently developed method in functional brain imag- 
ing analysis—connectome-based predictive modeling—to identify 
a brain network associated with high-creative ability, using func- 
tional magnetic resonance imaging (fMRI) data acquired from 
163 participants engaged in a classic divergent thinking task. At 
the behavioral level, we found a strong correlation between cre- 
ative thinking ability and self-reported creative behavior and ac- 
complishment in the arts and sciences (r — 0.54). At the neural 
level, we found a pattern of functional brain connectivity related 
to high-creative thinking ability consisting of frontal and parietal 
regions within default, salience, and executive brain systems. In a 
leave-one-out cross-validation analysis, we show that this neural 
model can reliably predict the creative quality of ideas generated by 
novel participants within the sample. Furthermore, in a series of ex- 
ternal validation analyses using data from two independent task fMRI 
samples and a large task-free resting-state fMRI sample, we demon- 
strate robust prediction of individual creative thinking ability from the 
same pattern of brain connectivity. The findings thus reveal a whole- 
brain network associated with high-creative ability comprised of cor- 
tical hubs within default, salience, and executive systems—intrinsic 
functional networks that tend to work in opposition—suggesting that 
highly creative people are characterized by the ability to simulta- 
neously engage these large-scale brain networks. 
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ehavioral and neuroimaging studies have begun to uncover 
the cognitive and neural processes that give rise to novel and 
useful ideas (1). Considerable research has focused on charac- 
terizing variation in creative thinking ability, often employing 
assessments of divergent thinking, which measure people’s ability 
to generate solutions to open-ended problems, such as inventing 
new uses for objects (2). Because divergent thinking ability 
moderately predicts real-world creative achievement in the arts 
and sciences (3), such laboratory-based assessments are thought 
to provide a reliable index of general creative ability. However, 
despite progress in the psychology and neuroscience of creativity, 
the field still lacks clarity on the neurocognitive characteristics 
that distinguish the highly creative brain (4). The present re- 
search thus aims to discover whether a specific brain connectivity 
profile characterizes high-creative thinking ability and to de- 
termine whether individual creativity can be reliably predicted 
from the strength of functional connectivity within this network. 
Using functional magnetic resonance imaging (fMRI), several 
studies have identified discrete brain regions that support per- 
formance on verbal creativity tasks, largely localized within 
frontoparietal and frontotemporal cortices involved in cognitive 
control and semantic memory retrieval, respectively (5, 6). More 
recently, researchers have embraced new techniques in {MRI 
data analysis to examine coordinated patterns of neural activity 
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across multiple distributed brain regions (i.e., functional con- 
nectivity) during various tasks that assess creative cognition 
and artistic performance, including divergent thinking, figurative 
language production, musical improvisation, poetry composition, 
and visual art production (1). This work highlights the contri- 
bution of three large-scale brain systems that dynamically in- 
teract to support creative task performance: the default mode 
network, comprised of cortical midline and posterior inferior 
parietal regions; the executive (or frontoparietal) control net- 
work, comprised of lateral prefrontal and anterior inferior pa- 
rietal regions; and the salience network, comprised of bilateral 
insula and anterior cingulate cortex (7). Building on dual-process 
theories of creative cognition (8), which emphasize idea gener- 
ation and evaluation processes, a recent brain network model 
posits that the default network supports idea generation and the 
executive network supports idea evaluation (1, 9), consistent with 
their established roles in mental simulation and executive cog- 
nition, respectively (7). Furthermore, the salience network— 
which contributes to the detection of behaviorally relevant 
stimuli and facilitates dynamic transitions between default and 
executive systems (10)—may identify candidate ideas stemming 
from generative processes within the default network and for- 
ward such information to executive systems for higher order 
processing (9). 


Significance 


People's capacity to generate creative ideas is central to tech- 
nological and cultural progress. Despite advances in the neuro- 
science of creativity, the field lacks clarity on whether a specific 
neural architecture distinguishes the highly creative brain. Using 
methods in network neuroscience, we modeled individual cre- 
ative thinking ability as a function of variation in whole-brain 
functional connectivity. We identified a brain network associ- 
ated with creative ability comprised of regions within default, 
salience, and executive systems—neural circuits that often work 
in opposition. Across four independent datasets, we show that a 
person's capacity to generate original ideas can be reliably 
predicted from the strength of functional connectivity within 
this network, indicating that creative thinking ability is charac- 
terized by a distinct brain connectivity profile. 


Author contributions: R.E.B., M.B., A.F., J.Q., T.R.K., M.J.K., and P.J.S. designed research; 
R.E.B., A.P.C., M.B., Q.C., A.F., J.Q., and P.J.S. performed research; M.D.R. contributed new 
reagents/analytic tools; R.E.B., Y.N.K., A.P.C., and M.D.R. analyzed data; and R.E.B., Y.N.K., 
A.P.C., M.D.R., M.B., Q.C., T.R.K., M.J.K., and P.J.S. wrote the paper. 


The authors declare no conflict of interest. 


This article is a PNAS Direct Submission. O.S. is a guest editor invited by the 
Editorial Board. 


Published under the PNAS license. 
‘To whom correspondence should be addressed. Email: rbeaty@fas.harvard.edu. 


This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. 
1073/pnas.1713532115/-/DCSupplemental. 


PNAS Early Edition | 1of6 


ul 
[e] 
z 
w 
9 
Nn 
e 
& 
> 
ul 
z 


PSYCHOLOGICAL AND 
COGNITIVE SCIENCES 


(cl NN 


L 


2N 


31 


A recent investigation of functional connectivity during di- 
vergent thinking (11) reported coupling between default and 
salience regions that preceded coupling between default and 
executive regions, potentially reflecting dynamic shifts in idea 
generation and evaluation across time. Moreover, participants 
that produced more original responses in a divergent thinking 
task showed higher global efficiency (i.e., a smaller number of 
paths needed to traverse between a given pair of brain regions) 
within a network comprised of default, salience, and executive 
network nodes. This pattern of functional connectivity has been 
observed across several task contexts requiring creative idea 
production, including studies with professional artists engaged in 
domain-specific tasks in the scanner (12). Functional coupling of 
these large-scale systems is particularly notable in light of past 
work reporting a negative or “anticorrelated” relation between 
the default and executive networks during task-free resting-state 
fMRI (13), as well as various task-based paradigms, including the 
well-documented finding of default network deactivation during 
executively demanding cognitive tasks (14). 

In the present research, we hypothesized that individual varia- 
tion in the ability to simultaneously engage the default, executive, 
and salience brain systems may provide a neurophysiological 
marker of creative thinking ability. We tested this hypothesis in a 
sample of 163 participants engaged in a creative thinking task 
during fMRI, and we employed connectome-based predictive 
modeling (cpm)—a method recently developed to predict aspects 
of human behavior (e.g., cognitive abilities) from patterns of 
whole-brain functional connectivity (15-17)—to examine whether 
creative thinking ability can be reliably predicted from an indi- 
vidual’s unique pattern of brain connectivity. The predictive power 
of cpm has recently been demonstrated in studies of fluid in- 
telligence (15) and sustained attention (18), revealing reliable 
prediction of these behavioral variables in participants whose data 
were not used in model construction (16). Here, we aimed to 
uncover a “creative connectome"—a whole-brain network asso- 
ciated with creative thinking ability—and conduct a leave-one-out 
cross-validation analysis (i.e., internal validation) to test whether 
the strength of functional connectivity within this network can 
reliably predict creative ability in novel participants within this 
dataset. To further assess the predictive power of this neural 
model, we conducted three external validation analyses using both 
task fMRI and resting-state fMRI data obtained from two dif- 
ferent laboratories and tested whether the strength of functional 
connectivity within this connectome can predict creative behavior 
in independent samples. 


Results 


Network Definition and Neuroanatomy. Functional imaging data 
were acquired as participants completed a classic divergent 
creative thinking task, where participants are presented with an 
object cue and are required to generate an unusual and creative 
use for it (Materials and Methods). The creative quality of re- 
sponses produced during the {MRI task, as well as two divergent 
thinking tasks completed outside the scanner, was assessed by 
four trained raters on a five-point scale (19). Using latent vari- 
able modeling (Fig. S1), we found that creativity scores strongly 
correlated with self-reported creative behavior and achievement 
in the arts and sciences [r(161) — 0.54, P « 0.001]—consistent 
with past work (20)—suggesting that these scores reflected an 
ecologically valid marker of creative ability. 

Whole-brain functional networks were constructed for each 
participant by extracting and correlating the task-related blood 
oxygen level-dependent (BOLD) signal from 268 brain regions 
(21). To identify network edges (ie., functional connections) 
related to creative thinking ability, we correlated all edges in this 
network with individual creativity scores extracted via latent 
variable modeling and applied a statistical threshold (P « 0.01) 
to retain the most significant edges in the connectivity matrices 
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(15, 16). This analysis revealed 224 edges that positively corre- 
lated with creative ability (“high-creative network") and 603 
edges that negatively correlated with creative ability ("low- 
creative network"; total possible edges — 35,778). Note that 
both low- and high-creative networks include some positive and 
negative connections, but the intersubject variability of these 
connections is either positively or negatively associated with 
behavior, respectively. 

The high-creative network exhibited dense functional con- 
nections in predominantly frontal and parietal cortices (Fig. 1). 
Consistent with past work, the regions showing the highest de- 
gree (k ; i.e., number of functional connections) corresponded to 
the core hubs of the default [e.g., left posterior cingulate cortex; 
Brodmann area (BA) 23; k = 19], salience (e.g., left anterior 
insula; BA 45; k = 14), and frontoparietal/executive network 
(e.g., right dorsolateral prefrontal cortex; BA 9; k = 8; Table S1 
and Dataset S1). Of the 25 highest degree nodes in the high- 
creative network, 12 were within the default network, 4 were 
within the salience/cingular-opercular network, and 3 were 
within the frontoparietal/executive network. The low-creative 
network showed diffuse connections across the brain, most 
prominently within subcortical/brainstem structures (e.g., left 
thalamus; k = 14), the default network (e.g., left posterior cingulate 
cortex; BA 23; k = 23), and cerebellum (Table S2 and Dataset S2). 
Of the 25 highest degree nodes in the low-creative network, 8 were 
within subcortical/brainstem structures, 5 were within the default 
network, and 4 were within the cerebellum. 


Internal Validation: Prediction of Creative Ability Using Task fMRI 
Data. We employed a leave-one-out cross-validation analysis 
(i.e., internal validation) to test whether the brain connectivity 
model (i.e., strength of functional connectivity within the high- 
and low-creative networks, respectively) (15 and 18) could re- 
liably predict creative thinking ability in novel participants (17). 
(Note that high- and low-creative networks can differ in each 
round of cross-validation since they are defined on a different set 
of 162 individuals and tested on the left-out 163rd participant.) 
Significant correlations emerged between the model-predicted 
and observed creativity scores in both the high-creative net- 
work [r(161) = 0.30, P = 1.29 x 1074] and low-creative network 
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Fig. 1. Depictions of the high- and low-creative networks. Circle plots (A) 
and glass brains (B) were thresholded to show the highest degree (k) nodes 
in the networks (high-creative k = 10, low-creative k = 18). Colors within the 
circle plots correspond to lobes of the brain. L, left hemisphere; R, right 
hemisphere. 
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[r(161) = 0.22, P = 0.004] (Fig. 2); a Steiger’s Z test for dif- 
ferences between two dependent correlations revealed that the 
correlations were not significantly different from one another 
(z = 0.83, P = 0.40) (22). 

To account for the nonindependence of the leave-one-out 
folds, we conducted a permutation test by randomly shuffling the 
creativity scores 5,000 times and rerunning the prediction pipe- 
line, to create a null distribution of r values for both high- and 
low-creative networks (18). The P values of the empirical cor- 
relation values, based on their corresponding null distribution, 
were computed by the following formula: (1 + the number of 
permutated r values greater than or equal to the empirical r)/ 
1,001. The permutation test results revealed significant P values 
for the empirical r values for both high-creative (P = 0.001) and 
low-creative (P = 0.019) networks. 

As a test of model specificity, we reran the prediction model 
using functional connectivity data from the noncreative control 
task [i.e., object characteristics task (OCT)] and the same crea- 
tive ability scores. Here, the model construction procedure was 
the same as above (i.e., trained on n — 1 participants’ creativity 
data), but the model was tested on the left-out participant's OCT 
connectivity data instead of their alternate uses task (AUT) 
connectivity data. The results revealed a weak, nonsignificant 
correlation between the observed and predicted creative-ability 
score for the high-creative network [r(161) = 0.05, P = 0.47], but 
a moderate and significant correlation between the observed and 
predicted creative-ability score for the low-creative network 
[r(161) = 0.19, P = 0.01]. Permutation tests using 1,000 iterations 
yielded a nonsignificant effect for the high-creative network (P — 
0.25) and a significant effect for the low-creative network (P — 
0.01), indicating that a model trained on the creative task data 
and tested on the control task data is moderately predictive of 
low-creative ability. 

In sum, the internal validation confirmed that individual dif- 
ferences in creative thinking ability can be reliably predicted in 
novel participants based on the strength of functional connec- 
tivity within task-related brain networks. Next, we further tested 
the predictive potential of the high- and low-creative network 
models by conducting three external validation analyses using 
fMRI data obtained from three independent samples of partic- 
ipants. Note that all three analyses testing associations with 
creative behavior use the same high- and low-creative network 
models derived above. 


External Validation 1: Prediction of Creative Ability Using Novel Task 
fMRI Data. We began by testing whether strength of connectivity 
within the high- and low-creative networks defined in the full set 
of 163 participants from the internal validation sample could 
predict creative thinking ability in a sample of Austrian partici- 
pants (n = 39) who performed a different divergent thinking task 
during {MRI (SI Materials and Methods). The high- and low- 
creative networks derived from the model training dataset above 
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Fig. 2. Correlation between predicted and observed creativity scores, 
standardized for visualization. 
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were applied to fMRI data acquired during an experiment that 
required participants to continuously generate creative object 
uses over an extended period (60 s). High- and low-creative 
network strength values were thus computed for each participant 
by summing correlation coefficients within the high- and low- 
creative networks. Results revealed a significant prediction of 
creativity scores (ie. composite ideational fluency and origi- 
nality) for the high-creative network [r (37) = 0.35, P = 0.03], but 
not the low-creative network [r (37) = —0.04, P = 0.78] (Fig. S2), 
indicating that participants with stronger functional connectivity 
within the high-creative network tended to generate more orig- 
inal ideas; a Steiger's Z test revealed that the correlations were 
not significantly different from one another (z = 1.64, P = 0.10). 
This external validation extends the internal validation by dem- 
onstrating the generalizability of the high-creative model to a 
novel and somewhat culturally distinct sample. 


External Validation 2: Prediction of Creative Ability Using Novel Task 
fMRI Data. A second external validation analysis was conducted 
with a new sample of Austrian participants who completed a 
divergent thinking task during fMRI (n = 54) (SI Materials and 
Methods). The high-creative and low-creative networks were 
again applied to each participant’s fMRI data to test whether 
creative thinking ability could be predicted from strength of 
connectivity within these networks. The strength of connectivity 
within the high-creative network was significantly correlated with 
behavioral performance (i.e., composite ideational fluency and 
originality) [r(52) = 0.28, P = 0.03] (Fig. S2). Performance was 
not significantly correlated with strength within the low-creative 
network [r(52) = 0.02, P = 0.91]. A Steiger's Z test revealed that 
the correlations were not significantly different from one another 
(z = 1.57, P = 0.11), consistent with the previous external 
validation analysis. 

To test whether the model selectively predicts performance on 
the cognitive task of interest (ie., divergent thinking) and not 
tasks that do not require divergent thinking, we ran a second 
analysis with functional connectivity matrices extracted from 
fMRI data during the control task in this experiment (i.e., gen- 
erating adjectives in response to word cues). The model was 
sensitive to the cognitive process engaged: Neither high-creative 
network strength [r(52) = 0.12, P = 0.38] nor low-creative net- 
work strength [r(52) = 0.03, P = 0.84] was significantly related to 
divergent thinking performance when based on functional con- 
nectivity data from the (noncreative) control task. Stated dif- 
ferently, the high- and low-creative networks applied to {MRI 
data during the control task were not significantly related to 
creative thinking ability. Likewise, when performance data from 
the noncreative control task was considered (i.e., mean num- 
ber of adjectives generated), neither high-creative network 
strength [r(52) = 0.11, P = 0.41] nor low-creative network 
strength [r(52) 2 0.13, P — 0.33] was related to performance: 
The high- and low-creative networks applied to fMRI data 
during the noncreative control task were not significantly re- 
lated to performance on the noncreative control task. To- 
gether, these findings provide evidence for the specificity of the 
high-creative network in predicting creative thinking ability. 


External Validation 3: Prediction of Creative Ability Using Resting- 
State fMRI Data. As an even more powerful test of external val- 
idity, we examined whether strength of connectivity within the 
high- and low-creative networks could predict creative thinking 
ability using task-free resting-state {MRI data. Given past work 
showing a correspondence between resting-state and task-in- 
duced network architecture (23), we hypothesized that strength 
of functional connectivity within the high- and low-creative 
networks at rest would predict creative ability, albeit to a lesser 
extent than with task-state connectivity data given the model's 
apparent sensitivity to task-congruent neural activity. 
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We obtained resting-state fMRI data from a large sample of 
Chinese participants (n — 405) who completed a battery of di- 
vergent thinking tasks outside of the scanner (SI Materials and 
Methods). The high-creative and low-creative networks were 
again applied to the resting-state fMRI data of the participants. 
Creativity scores correlated significantly with high-creative net- 
work strength [r(403) = 0.13, P = 0.008], as well as low-creative 
network strength [r(403) = 0.11, P = 0.03] (Fig. S2), indicating 
that the model can reliably predict creative performance from 
intrinsic functional brain connectivity. (Note that the positive 
relationship between creative performance and low-creativity 
network strength is unexpected. Given the failure of the low- 
creativity network to generalize to external validation sets 1 and 
2, this network may be a less reliable predictor of creativity than 
the high-creativity network.) A Steiger's Z test revealed that the 
correlations were not significantly different from one another 
(z = 0.79, P = 0.42). As a test of model sensitivity, we assessed 
whether high- and low-creative network strength could predict 
performance on a test of fluid intelligence. Neither high-creative 
[r(403) = 0.04, P = 0.39] nor low-creative [r(403) = 0.04, P = 
0.40)] network strength correlated with intelligence scores, thus 
demonstrating that the neural models selectively predict creative 
thinking ability and not simply higher cognitive ability. 


Discussion 


Using a recently developed technique in fMRI data analysis— 
connectome-based predictive modeling (15-17)—we uncovered 
whole-brain networks associated with performance on a classic 
assessment of general creative thinking ability. This high-crea- 
tive network exhibits dense functional connections between 
core nodes of the default, executive, and salience systems— 
networks that typically work in opposition—suggesting that the 
creative brain is marked by a tendency to simultaneously en- 
gage these large-scale circuits to a greater degree than the less 
creative brain. Critically, we demonstrate the robustness of the 
model in predicting creative ability across three independent 
datasets obtained from two different laboratories, including two 
task-based {MRI samples and one task-free resting-state sam- 
ple. Taken together, the findings provide evidence that creative 
thinking ability can be reliably predicted from an individual's 
unique brain connectivity profile. 

The neuroanatomy of the high-creative network is remarkably 
similar to patterns of functional connectivity reported in recent 
studies of creative cognition (1). Specifically, we found that the 
regions showing the greatest number of significantly correlated 
functional connections corresponded to the hubs of three large- 
scale brain networks: default (posterior cingulate cortex), exec- 
utive (right dorsolateral prefrontal cortex), and salience (left 
anterior insula). These regions showed increased functional 
connectivity in a recent {MRI study of divergent thinking (11), 
which also reported a positive correlation between individual 
creativity scores and global efficiency within a set of regions 
consisting of the core nodes of these three networks, consistent 
with the connectivity profile associated with high-creative ability 
in the present study. Notably, neuroimaging studies of artistic 
performance have also reported functional coupling among re- 
gions within these networks (12, 24, 25), suggesting that the 
ability to coactivate brain systems that tend to work in opposition 
at rest and during demanding cognitive tasks (14) may reflect a 
domain-general mechanism of creative information processing. 

The low-creative network was largely comprised of subcortical, 
cerebellar, and sensorimotor brain regions. Although we did not 
have a priori hypotheses about the low-creative network, a few 
points are worth noting. Regarding neuroanatomy, the regions 
showing the most connections within the network are less com- 
monly reported in neuroimaging studies of creative cognition (5, 
6), suggesting that they are less relevant to successful perfor- 
mance on creative thinking tasks. One possibility, then, is that 
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low-creative ability is characterized by increased interactions 
among brain regions that do not reliably support creative cogni- 
tion. Although the network included some regions of the default 
and salience systems, a closer examination of the connectivity 
pattern showed that these systems were not functionally linked— 
as in the high-creative network—but rather showed greater 
coupling with subcortical, sensorimotor, and cerebellar regions. 
Previous evidence has found that these regions support procedural 
(habitual) responses, which suggests that low-creative people 
might be retrieving cached or previously learned responses (26). 
In light of past work linking default activity to strong semantic 
associations (27) and automatic response tendencies (28), low- 
creative ability may be characterized by an inability to transcend 
salient conceptual knowledge when drawing on memory repre- 
sentations within the default network. Thus, when attempting to 
generate original ideas, low-creative individuals may activate 
common, cached semantic associations that are, in turn, not ef- 
fectively regulated by the executive network, which may function 
to inhibit prepotent response tendencies and redirect search 
processes in the highly creative brain (29). 

Moreover, although the low-creative network showed signifi- 
cant performance prediction in the original dataset, the external 
validation analyses only revealed significant correlations with 
creative thinking ability in one of the three samples (i.e., external 
validation 3; resting-state {MRI)—and this relationship was in 
the unexpected direction—suggesting that the network may be 
less reliable in its prediction of low-creative ability. It is impor- 
tant to note, however, that, across all datasets, the numerically 
different correlations between creativity scores and high- versus 
low-creative network strength were not significantly different 
from one another (and our samples were underpowered to de- 
tect such differences). Nevertheless, the high-creative network 
showed stable prediction of creative ability across datasets, and 
the neuroanatomy of this network mirrors previous fMRI studies 
of creative cognition reporting functional interactions among 
default, executive, and salience networks (1), suggesting that it 
may be a more reliable neural marker of creative ability in future 
work. It is also worth noting that, although we found stronger 
prediction of creative thinking ability using task-based fMRI 
data, this pattern also emerged using resting-state data, in- 
dicating that the brain-behavior correlations are driven in part by 
stable trait-level variation in functional connectivity. Thus, both 
intrinsic network structure as well as task-induced network re- 
organization likely capture important individual variation (30). 

Coordination of the default, executive, and salience networks 
is consistent with a recent model of brain dynamics supporting 
creative cognition (1). According to this framework, the default 
network contributes to the generation of ideas via flexible and 
spontaneous combinatory mechanisms involved in memory re- 
trieval and mental simulation. The salience network, in turn, 
functions to identify candidate ideas— potentially useful in- 
formation generated via the default network—and forward such 
information to frontoparietal executive systems for high-order 
processing (e.g., idea evaluation, elaboration, or revision). Al- 
though the model implies a serial progression of idea generation 
and evaluation, recent evidence suggests that executive systems 
may interact with ongoing generative processes within the de- 
fault network by imposing constraints on performance and 
maintaining higher order goals (1). Goal maintenance may thus 
benefit creative performance by guiding and constraining spon- 
taneous cognition to meet specific creative goals (31). 

Creativity remains a complex construct that will require con- 
siderable further research to uncover its many manifestations in 
the brain. Unlike some aspects of cognition that have been re- 
liably localized to specific brain regions, complex constructs such 
as creativity are likely a product of similarly complex neural 
mechanisms that engage the whole brain. Thus, in addition to 
pursuing classic questions related to where creativity occurs in 
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the brain, a promising direction for future research will be to focus 
on how the brain thinks creatively. Network-based approaches are 
particularly well-suited to address such questions because they can 
accommodate the complex interplay of multiple neurocognitive 
processes (e.g., memory retrieval, mental simulation, and cogni- 
tive control). We encourage future research to further explore 
how the brain thinks creatively across different task contexts and 
to develop behavioral and neural interventions (e.g., using brain 
stimulation) to enhance creative performance by targeting large- 
scale networks and their dynamic interactions. 


Materials and Methods 


Participants. The total sample consisted of 163 participants recruited from the 
University of North Carolina at Greensboro (UNCG) and the surrounding 
community (113 women, mean age = 22.50 y, SD = 5.79) and specifically 
over-sampled art, music, and science majors to increase the sample's pop- 
ulation of creative domains. Participants were recruited as part of a larger 
study on individual differences in creativity (which involved numerous lab- 
oratory and ecological measures and procedures not discussed here) and 
were paid up to $100 for their time. All participants were right-handed with 
normal or corrected-to-normal vision and reported no history of neurolog- 
ical disorder, cognitive disability, or medication that affects the central 
nervous system. Participants provided written informed consent. The study 
was approved by the UNCG Institutional Review Board. 


fMRI Task Paradigm. Participants completed a creativity task and a control 
task in an event-related design during functional imaging: an alternate 
uses task (AUT) of divergent thinking and an object characteristics task (OCT). 
The task paradigm was similar to protocols used in past research (11, 32). In 
the AUT, participants were presented with a common object cue and asked 
to imagine a new and unusual use for it; in the OCT, participants were 
presented an object cue and asked to think of typical object characteristics. 
Participants received thorough training on both tasks and completed several 
practice trials before scanning. 

The task paradigm consisted of a jittered fixation cross (4-6 s), a cue in- 
dicating the upcoming condition ("create" or "object"; 3 s), a thinking period 
presenting an object cue in text (12 s), and a verbal response period requiring 
participants to verbalize their response into an MRI-compatible microphone 
(5 s). In the AUT, participants were asked to use the thinking period to 
imagine creative uses for the object and then speak their most creative re- 
sponse into the microphone during the response period (29, 33); in the OCT, 
participants were asked to think of the most common physical characteristic 
of the object. Consistent with past work, an experimenter recorded verbal 
responses in real time for subsequent analysis of idea quality. 


Behavioral Assessment. In-scanner performance on the AUT was assessed by 
recording verbal responses via MRI-compatible microphone, which were 
subsequently coded for creative quality. Participants also completed two 
AUTs as part of the postscan behavioral assessment to obtain data on their 
performance in a standard testing environment. Consistent with conven- 
tional procedures, participants were given 2 min to continuously generate 
alternate uses for each of two objects (i.e., box and rope), which were 
subsequently coded for creative quality by four trained raters along with the 
in-scanner responses, using a 1 (not at all creative) to 5 (very creative) scale; 
raters were instructed to consider uncommonness, remoteness, and clever- 
ness when coding responses, but to provide a single holistic score for each (11, 
19, 34, 35). The raters scored all in-scanner responses (23 trials) and labora- 
tory-based responses (two tasks). For in-scanner responses, the average of 
the scores for each rater was computed, resulting in good interrater re- 
liability [intraclass correlation (ICC) = 0.78]; for laboratory-based responses, 
the ratings for each task were similarly averaged, also resulting in good rater 
agreement in the box (ICC = 0.69) and rope tasks (ICC = 0.75). Finally, a 
battery of questionnaires was administered to assess creative achievement 
and everyday creative behavior: the Creative Achievement Questionnaire 
(CAQ) (36); the Biographical Inventory of Creative Behavior (BICB) (37); and 
the Inventory of Creative Activities and Achievements (ICAA) (38). 

Latent variable modeling was used to model the creativity ratings of the in- 
scanner and laboratory-based AUT responses. We specified a higher order 
latent variable model using Mplus 7.2. The higher order creativity factor was 
indicated by two lower order latent factors: in-scanner ratings and labora- 
tory-based ratings. The in-scanner factor was indicated by four observed 
variables (i.e., the average creativity ratings of the four raters). The labo- 
ratory-based factor was indicated by two lower-order factors representing 
the two laboratory-based tasks, which were both in turn indicated by four 
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observed variables (i.e., the average creativity ratings of the four raters; see 
Fig. S1). A “creative behavior and achievement" variable was also modeled, 
indicated by four observed variables, to assess the latent correlation between 
creative thinking ability and real-world creative behavior and accomplish- 
ment. For model identification, the paths of all indicators were constrained to 
equality, and the variances of the latent variables were fixed to 1. 


MRI Data Acquisition and Preprocessing. Participants completed the tasks in a 
single fMRI run. Whole-brain imaging was performed on a 3T Siemens 
Magnetom MRI system (Siemens Medical Systems) using a 16-channel head 
coil. BOLD-sensitive T2*-weighted functional images were acquired using a 
single shot gradient-echo echo-planar imaging (EPI) pulse sequence [repe- 
tition time (TR) = 2,000 ms, echo time (TE) = 30 ms, flip angle = 78°, 32 axial 
slices, 3.5 x 3.5 x 4.0 mm, distance factor 0%, field of view (FoV) = 192 x 
192 mm, interleaved slice ordering] and corrected online for head motion. 
The first two volumes were discarded to allow for T1 equilibration effects. 
Visual stimuli were presented using E-Prime and viewed through a mirror 
attached to the head coil. In addition to functional imaging, a high resolu- 
tion T1 scan was acquired for anatomic normalization. Functional volumes 
were slice time-corrected and realigned using the Statistical Parametric 
Mapping (SPM) 12 package (six motion parameters), coregistered and nor- 
malized to the Montreal Neurological Institute (MNI) template brain, and 
smoothed with an 8-mm? isotropic Gaussian kernel. To determine whether 
creativity scores (i.e., latent variable factor scores) were associated with 
participant head motion during functional imaging, we correlated these 
scores with mean frame-to-frame displacement. Results showed that crea- 
tivity scores were marginally but not significantly related to head motion 
(r = —0.14, P = 0.07). 

Task-related functional connectivity was assessed using the CONN toolbox 
(https://www.nitrc.org/projects/conn) (39) in MATLAB. For each participant, 
CONN implemented CompCor, a method for identifying principal compo- 
nents associated with segmented white matter (WM) and cerebrospinal fluid 
(CSF). In a first-level analysis, CompCor components and first-order deriva- 
tives of motion were entered as confounds and regressed from the BOLD 
signal. Additional preprocessing steps included high-pass filtering, linear 
detrending, and regression of outlying functional volumes (>97th percentile 
in normative sample; global-signal z-value threshold = 5, subject-motion mm 
threshold = 0.09) identified using the artifact removal toolbox (ART) (https:// 
www.hitrc.org/projects/artifact detect/) We also regressed the task struc- 
ture corresponding to the onsets and durations of the verbal response pe- 
riods to account for expected artifacts related to participant vocalization. 
Because CompCor can account for subject movement effects and other 
sources of noise in the BOLD signal, the global signal was not regressed. 


Functional Network Construction. Whole-brain networks were computed for 
each participant using CONN. Consistent with past work employing cpm (15, 
18), we used the Shen brain atlas, which consists of 268 regions of interest 
(ROIs) of 2-mm dimensions and provides whole-brain coverage of the cere- 
bral cortex, cerebellum, and brainstem (21). BOLD signal was extracted from 
each ROI during the thinking period of the AUT (23 trials, 12 s; collapsing 
across trials), and bivariate correlations were computed between each pair 
of ROls, resulting in a 268 x 268 correlation matrix for each participant. Note 
that the same anatomical parcellation and procedure were used to construct 
functional networks for all external validation analyses described here. 


Connectome-Based Predictive Modeling. The main analysis employed cpm to 
estimate participants' creative thinking ability from whole-brain, task- 
related functional connectivity. Cpm is a recently developed method for 
identifying functional brain connections related to a behavior variable of 
interest, which are then used to predict behavior in novel participants (i.e., 
participants whose data were not used in model creation) (17). The cpm 
procedure was recently described in a series of studies reporting its appli- 
cation to cognitive variables such as fluid intelligence and attention control 
(18), but we briefly summarize the cpm processing pipeline here. The 
MATLAB syntax used for cpm is freely available online (https://www.nitrc. 
org/projects/bioimagesuite/), as is the connectome visualization software 
(bisweb.yale.edu/build/connviewer.html). 

In a first step, a vector of behavioral values (i.e., a single creativity factor 
score for each participant) was correlated with each edge (i.e., correlation of 
mean BOLD signals between a given pair of brain regions) in the functional 
connectivity matrix of each participant. Next, a threshold was applied to the 
matrix to retain only edges that were significantly positively and negatively 
correlated with behavior (P « 0.01). Cpm was then applied to a participant's 
data by summing the edge strength (i.e., correlation coefficients) in the 
positive and negative tails of correlation distribution; the frequency 
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distributions of behavioral and connectivity values were checked for nor- 
mality to meet assumptions for Pearson correlations. Next, a linear regression 
model was specified to estimate the relationship between the model predicted 
behavior score and the observed behavior score. Finally, the model was applied 
to new participants in a level-one-out cross-validation, such that the model was 
trained on n — 1 participants’ connectivity matrices and behavior scores, and 
tested on the left-out participant. Note that feature selection (i.e., network 
edges retained) also occurs within the leave-one-out loop, resulting in slightly 
different networks and predictive models for each iteration. The predictive 
power of the resulting model is reflected in the magnitude and statistical sig- 
nificance of the Pearson correlation between the model predicted and ob- 
served behavior scores. For a thorough description of cpm, see Shen et al. (17). 


External Validation. One hundred sixty-three different positive and negative 
networks and linear models—one for each round of leave-one-out cross- 
validation—were defined during internal validation. To generate a final 
model to apply to completely independent samples, we defined final high- 
and low- creativity networks using data from all 163 training subjects. The 
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"final" high creativity network included 224 edges, and the final low crea- 
tivity network included 603 edges. These networks are visualized in Fig. 1. 
The networks were highly similar to networks identified in leave-one-out 
folds during internal validation: On average, edges in the final high-creative 
network appeared in 151 of the 163 positive masks generated during leave- 
one-out cross-validation, and edges in the final low-creative network 
appeared in 154 of the 163 negative masks. Next, we calculated strength in 
these networks for each subject in three completely independent datasets. 
We correlated network strength and behavior in these novel samples, re- 
vealing that, in each case, individuals with higher high-creative networks 
were more creative. 
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SI Results 


Background Connectivity Analyses. Previous research has indicated 
that coactivations of brain regions can increase correlations (i.e., 
functional connectivity) by introducing a common timing struc- 
ture (1). We thus sought to verify that our model does not rely on 
task-evoked coactivation patterns to predict individual differ- 
ences in creativity by performing analyses on background con- 
nectivity. Specifically, we modeled the effect of task using the 
CONN toolbox and calculated functional connectivity matrices 
using the residuals of this regression. All other preprocessing 
steps were identical to those described in the main text. 
Performing leave-one-subject-out cross-validation in our initial 
training data set with these new connectivity matrices, we found 
similar predictive power for the high-creativity network [r(161) — 
0.32, P = 2.76e—05] but not the low-creativity network [r(161) = 
0.15, P = 0.05], which in the main text we describe as less reli- 
able. The networks themselves showed considerable overlap with 
those visualized in Fig. 1. External validation revealed that the 
high-creativity network model generalized to predict creativity 
from functional connectivity in both the first external task vali- 
dation sample [n = 39; r(37) = 0.47, P = 0.002], as well as the 
second [n = 54; r(52) = 0.32, P = 0.02], and generalized to predict 
creativity in a third independent group from resting-state func- 
tional connectivity [n = 405; r(403) = 0.10, P = 0.04]. Numeric 
differences in the results do not affect our main conclusion that 
models based on patterns of functional connectivity predict indi- 
vidual differences in creativity in novel individuals and datasets. 


Partial Correlation Analyses. The Pearson correlations reported in 
the main text reflect associations between creative thinking ability 
and network connectivity strength. However, because creativity 
scores were moderately correlated with in-scanner motion in the 
original dataset, we sought to determine whether brain-behavior 
relations reported in the text were robust to subject motion. Here, 
we hence report partial correlations between observed and 
predicted creativity scores in the original dataset, and correlations 
between creativity scores and network strength in the external 
validation datasets—controlling for mean in-scanner motion 
values. Regarding the original dataset, we found that controlling 
for motion resulted in slightly attenuated correlations between 
observed and predicted creativity scores in the high-creative 
[r(160) = 0.28, P < 0.001] and low-creative [r(160) = 0.19, P = 
0.01] networks. In the first task validation dataset (n — 39), we 
found a slight increase in the correlation between creativity 
scores and network strength in the high-creative network [r(36) — 
0.42, P = 0.009] and a similarly small and nonsignificant corre- 
lation between creativity scores and low-creative network 
strength [7(36) = —0.08, P = 0.65]. Regarding the second task 
validation dataset (n — 54), we found a comparable pattern of 
correlations as previously reported for both the high-creative 
network [r(51) = 0.29, P = 0.04] and low-creative network 
[r(51) = —0.02, P = 0.91]. Likewise, regarding the resting-state 
dataset (n — 405), the previously reported effects were largely 
similar for both the high-creative network [r(402) = 0.13, P = 
0.01] and low-creative network [r(402) = 0.10, P = 0.04]. 


SI Materials and Methods 


External Validation 1: Task fMRI. 

Participants. Behavioral and functional imaging data were rean- 
alyzed from a previous study that assessed brain activity during the 
AUT (2). The total sample included 39 participants (26 female, 
mean age = 22.50 y, SD = 3.12). All participants were right- 
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handed with normal or corrected-to-normal vision and reported 
no history of neurological disorder, cognitive disability, or medi- 
cation that affects the central nervous system. Participants com- 
pleted consent forms and were paid for their time. The study was 
approved by the local ethics committee of the University of Graz. 
fMRI task paradigm. Participants completed the AUT during func- 
tional imaging. Similar to the above design and previous studies, 
the task involved generating unusual and creative uses for everyday 
objects. A unique feature of this paradigm, however, was that it 
allowed continuous idea generation over extended trial durations 
(i.e., 1 min vs. 12s in study 1): Participants spontaneously generated 
ideas in a self-paced fashion and vocalized them as soon as they 
came to mind. Functional imaging data were acquired in a single 
run that included 15 task blocks. A trial block consisted of a jittered 
fixation cross (20—22 s) followed by a thinking period presetting an 
object cue (60 s) where participants were required to think of 
creative uses. Participants continuously generated ideas during the 
thinking period and vocalized their responses into a mouthpiece 
fixed to the head coil; the mouthpiece led to a microphone lo- 
cated outside the scanning room where an experimenter tran- 
scribed the responses. After the scan, participants were presented 
with their responses and asked to label each as *old" (i.e., recalled 
from memory) or “new” (i.e., generated on the spot). In lieu of a 
control task, the study contrasted neural activity associated with 
the generation of old vs. new ideas. However, for the current 
analysis, we analyzed data from the entire functional run, con- 
sistent with the seminal studies using cpm (3). 

Behavioral assessment. Participants’ verbal responses to the AUT 
were scored for creative quality by three trained raters using a 
subjective scoring method; raters used a four-point scale to code 
responses (1 = uncreative, 4 = very creative) and were instructed 
to consider both the originality and appropriateness of each idea 
into a single rating. Analysis of rater agreement showed good 
interrater reliability (ICC — 0.78). A fluency score was also 
computed for each participant by summing the total number of 
ideas generated for each trial and computing an average across 
the 12 trials. In addition to the in-scanner task responses, par- 
ticipants also completed laboratory-based AUT measures. Par- 
ticipants were given 2 min to type unusual and creative uses for 
everyday objects. Responses were similarly scored for fluency 
and originality by trained raters. A composite AUT score was 
computed for each participant by averaging the in-scanner and 
laboratory-based originality and fluency ratings. 

MRI data acquisition and preprocessing. Participants completed the 
tasks in a single {MRI run. Whole-brain imaging was performed 
on a 3T Tim Trio system (Siemens Medical Systems) using a 32- 
channel head coil BOLD-sensitive T2*-weighted functional 
images were acquired using a single shot gradient-echo EPI 
pulse sequence (TR = 2,000 ms, TE = 25 ms, flip angle = 90°, 
34 axial slices, slice thickness = 3 mm, FoV = 192 x 192 mm, 
interleaved slice ordering) and corrected online for head motion. 
The first two volumes were discarded to allow for T1 equili- 
bration effects. A high resolution T1 scan was acquired for an- 
atomic normalization. Visual stimuli were presented using 
Presentation software (Neurobehavioral Systems) and viewed 
through a mirror attached to the head coil. Functional volumes 
were slice time-corrected and realigned using SPM 12, anatom- 
ically coregistered, normalized to the MNI template brain 
(Montreal Neurological Institute), and smoothed with a 6-mm? 
isotropic Gaussian kernel. [Note that smoothing kernels var- 
ied across datasets because, in some cases, only preprocessed 
data (i.e., smoothed images) were available for reanalysis.] 


1 of 6 


Task-related functional connectivity was assessed using the 
CONN toolbox (4) in MATLAB using the same preprocessing 
procedures described above, including regression of verbal re- 
sponse onsets and durations. For functional network construction, 
the mean BOLD signal was extracted from each ROI during the 
thinking period of the AUT (12 trials, 60 s; collapsing across tri- 
als). [Whole-brain correlation matrices were computed for all 
external validation analyses presented here, but only those nodes 
and edges in the high- and low-creative networks (ie., those 
retained after statistical thresholding) were considered when 
computing network strength values in the validation samples.] 


External Validation 2: Task fMRI. 

Participants. Behavioral and functional imaging data were rean- 
alyzed from a previous study that assessed brain activity during the 
AUT (5). The sample included 54 participants (26 female, mean 
age = 24.06 y, SD = 2.99). All participants were right-handed with 
normal or corrected-to-normal vision and reported no history of 
neurological disorder, cognitive disability, or medication that af- 
fects the central nervous system. Participants completed consent 
forms and were paid for their time. The study was approved by the 
local ethics committee of the Medical University of Graz. 

fMRI task paradigm. Participants completed a creative thinking and 
control task during functional imaging: an AUT and an instances 
task (IT). In the AUT, participants were presented with a 
common object and required to generate unusual and original 
uses for it; in the IT—the control task—participants were pre- 
sented with an adjective and required to generate conditions or 
facts related to it (6); the AUT and IT trials were presented in a 
random order. In the current analysis, we were primarily in- 
terested in extracting neural activity during the AUT, but we also 
assessed activity during the IT to test the sensitivity of the pre- 
diction model to the cognitive process engaged. Participants 
continuously generated responses to the AUT or the IT during 
an idea generation period. Functional imaging data were ac- 
quired in a single run that included 40 trials of each condition. A 
trial consisted of a jittered fixation cross (4-8 s), a thinking pe- 
riod presenting a cue word (a noun in the AUT or an adjective in 
the IT; 15 s), and a verbal response period requiring participants 
to speak their responses into an MRI-compatible microphone 
(7 s); an experimenter recorded verbal responses online for 
subsequent analysis of idea quality. 

Behavioral assessment. Participants’ verbal responses to the AUT 
were scored for originality by four trained raters using a subjective 
scoring method; raters used a four-point scale to code responses 
(0 = not at all original, 3 = very original) and were instructed to 
consider both the originality and appropriateness of each idea into 
a single rating. Analysis of rater agreement showed good inter- 
rater reliability (ICC = 0.75). A fluency score was also computed 
for each participant by summing the total number of ideas gen- 
erated for each trial and computing an average across the 12 trials. 
A composite AUT score was computed for each participant by 
averaging the originality and fluency ratings. 

MRI data acquisition and preprocessing. Participants completed the 
tasks in a single {MRI run. Whole-brain imaging was performed on 
a 3T Siemens Magnetom Skyra (Siemens Medical Systems) using a 
32-channel head coil. BOLD-sensitive T2*-weighted functional 
images were acquired using a single shot gradient-echo EPI pulse 
sequence (TR = 2,400 ms, TE = 30 ms, flip angle = 90°, 36 axial 
slices, slice thickness = 3.5 mm, FoV = 240 mm) and corrected 
online for head motion. The first two volumes were discarded to 
allow for T1 equilibration effects. A high resolution T1 scan was 
acquired for anatomic normalization. Visual stimuli were pre- 
sented using Presentation software (Neurobehavioral Systems) 
and viewed through a mirror attached to the head coil. Functional 
volumes were slice time-corrected and realigned using SPM 8, 
normalized to the MNI template brain (Montreal Neurological 
Institute), and smoothed with an 8-mm? isotropic Gaussian ker- 
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nel; verbal response onsets and durations were also regressed 
from the BOLD signal. Functional connectivity was similarly 
assessed using the CONN toolbox (4) in MATLAB. As described 
above, the Shen anatomical atlas was used to define whole-brain 
networks. Here, we extracted mean BOLD signal from the 268 
ROIs during the AUT and IT separately (20 trials each, 15 s; 
collapsing across trials); bivariate correlations were computed 
between each pair of ROIs, resulting in a 268 x 268 correlation 
matrix for each of the 54 participants. 


External Validation 3: Resting-State fMRI. 

Participants. Resting-state functional imaging data were obtained 
from the Southwest University Longitudinal Imaging Multimodal 
(SLIM) brain data repository, a publicly available database of 
neuroimaging and phenotypic data acquired from a large sample 
of Chinese young adults (fcon 1000.projects.nitrc.org/indi/retro/ 
southwestuni qiu index.html) (7). The current sample included 
a large subset of participants who completed a battery of di- 
vergent thinking assessments, including the AUT (n = 405; 
215 females, mean age = 20.03 y, age range = 1.25). All par- 
ticipants were right-handed with normal or corrected-to-normal 
vision and reported no history of neurological disorder, cognitive 
disability, or medication that affects the central nervous system. 
Participants completed consent forms and were paid for their 
time. The study was approved by the Southwest University Brain 
Imaging Center Institutional Review Board. 

Behavioral assessment. The verbal form of the Torrance Tests of 
Creative Thinking (TTCT) (8) was used to assess creativity (i.e., 
divergent thinking ability). Five tasks from the Chinese version of 
the TTCT were administered: generating questions; causes and 
consequences; improving products; alternate uses; and manipu- 
lating objects. The TTCT provides a total creativity score as well 
as indices and scores for evaluating different creative processes or 
dimensions, which includes (i) fluency (the number of mean- 
ingful and relevant responses), (i) flexibility (the number of 
different categories of responses), and (iii) originality (the de- 
gree of originality of the responses). Three trained raters coded 
the creative quality of responses consistent with conventional 
scoring procedures described in the TTCT manual (ICC » 0.90). 
Latent variable analysis was applied to extract a latent variable 
creativity score for each participant. We specified the latent 
variable using Mplus 7.2, indicated by the average fluency, 
flexibility, and originality scores. To assess fluid intelligence, 
participants completed the Combined Raven's Test (CRT) 
(72 items), a widely used intelligence assessment in China that is 
based on Raven's Color Progressive Matrices. 

MRI data acquisition and preprocessing. Resting-state {MRI data were 
acquired for 8 min. Whole-brain imaging was performed on a 3T 
Siemens Trio MRI system (Siemens Medical Systems) using a 
12-channel head coil. BOLD-sensitive T2*-weighted functional 
images were acquired using a single shot gradient-echo EPI 
pulse sequence (TR = 2,000 ms, TE = 30 ms, flip angle = 90°, 
32 axial slices, 3.4 x 3.4 x 4.0 mm, FoV = 220 x 220 mm, in- 
terleaved slice ordering, 242 volumes) and corrected online for 
head motion. During functional imaging, participants were asked 
to keep their eyes closed, remain awake, and not think about 
anything in particular. A high resolution T1 scan was also ac- 
quired for anatomic normalization. Imaging data were pre- 
processed using SPM 12. The first 10 volumes from each 
subject's functional imaging data were discarded to account for 
steady-state magnetization, resulting in 232 volumes for sub- 
sequent analysis. Functional volumes were slice time-corrected, 
realigned, resliced to a voxel size of 3 mm?, normalized to the 
MNI template brain, and smoothed with a 6-mm? isotropic 
Gaussian kernel. Additional preprocessing included regression 
of motion parameters, along with white matter and cerebrospinal 
fluid (CSF); BOLD time series were filtered using a standard 
low-pass filter (0.01-0.08 HZ). 
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Fig. S1. Latent variable model of creative thinking and behavior. bicb, Biographical Inventory of Creative Behavior; caq, Creative Achievement Questionnaire; 
crea beh, creative behavior; crea cog, creative cognition (divergent thinking creativity); dt1, divergent thinking task 1 (laboratory responses); dt1 r1-r4, di- 
vergent thinking task 1 rater 1-4 scores (laboratory responses); dt2, divergent thinking task 2 (laboratory responses); ic act = Inventory of Creative Activities; 
ica ach = Inventory of Creative Accomplishments; mri r1-r4 = divergent thinking 1—4 scores (in-scanner responses). 
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Fig. S2. Scatterplots for all external validation analyses. High- and low-creative network strength values and observed creativity scores are standardized for 


visualization. Correlations between high- and low-creative network strength in external validation 1 (A; task fMRI), external validation 2 (B; task fMRI), and 
external validation 3 (C; resting-state fMRI). 
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Table S1. High-creative network neuroanatomy (top 25 high-degree nodes) 

No. K Node name Network L/R Lobe BA MNI 

1 19 Posterior cingulate Default L Limbic 23 —5, —36, 32 

2 15 Anterior insula Cingular-opercular L Insula 13 38.7, 8.1, —4.8 
3 14 Anterior insula Salience L Prefrontal 45 —32.5, 22.1, 5.8 
4 12 Precuneus Default L Limbic 31 —6.5, —53.9, 37.4 
5 9 Precuneus Default R Limbic 31 6.2, —57.4, 38.2 
6 8 Dorsolateral PFC Frontal-parietal R Prefrontal 10 44.6, 46.2, —4.9 
7 8 Anterior cingulate Default L Limbic 32 —6, 34.1, 26.3 
8 8 Posterior cingulate Default R Subcortical n/a 21.2, —36.4, 22.6 
9 8 Anterior temporal Default R Temporal 22 56.5, —8.5, —14.3 
10 7 Precentral gyrus Default L Prefrontal 8 —39.3, 1.7, 46.7 
11 7 Dorsolateral PFC Frontal-parietal L Prefrontal 9 46.1, 28.2, 26.8 
12 7 Ventromedial PFC Default L Prefrontal 11 —8.2, 39.7, -21.4 
13 7 Midcingulate n/a R Limbic 23 7, —18.8, 29.8 
14 6 Caudate n/a L Subcortical n/a —10.8, 23.9, 9.6 
15 6 Retrosplenial Default L Limbic 23 —8.6, —58.8, 17.6 
16 6 Inferior frontal gyrus Frontal-parietal L Prefrontal 44 —53.1, 18.4, 10.6 
17 6 Calcarine Visual R Limbic 23 28.4, —53.8, 7.1 
18 6 Posterior cingulate Default R Limbic 23 5.1, —38.9, 27 
19 6 Ventromedial PFC Default R Prefrontal 10 8.2, 45.9, —1.7 
20 5 Crus 1 Cerebellum L Cerebellum n/a —40.3, —74.2, —29.1 
21 5 Inferior frontal gyrus n/a L Prefrontal 47 —32, 20.5, —16 
22 5 Inferior frontal gyrus n/a L Prefrontal 11 —18.2, 19, -21 
23 5 Precuneus Salience R Limbic 31 8.3, —39.9, 48.1 
24 5 Precuneus Default R Limbic 23 12.3, —57.2, 18.1 
25 5 Anterior insula Salience R Insula 13 37.5, 21.1, —10.1 


BA, Brodmann area; K, degree; L, left; n/a, not available; PFC, prefrontal cortex; R, right. 


Table S2. Low-creative network neuroanatomy (top 25 high-degree nodes) 

No. K Region Network UR Lobe BA MNI 

1 23 Brainstem n/a L Brainstem n/a —7.2, —33, —39.4 
2 23 Posterior cingulate Default L Limbic 23 —5, —36, 32 

3 21 Precuneus Default L Limbic 31 —6.5, —53.9, 37.4 
4 21 Brainstem n/a R Brainstem n/a 7.5, —34.2, —37.3 
5 19 Fusiform gyrus Visual L Occipital 19 —43.2, —70.4, —13.8 
6 19 Brainstem n/a R Brainstem n/a 6, —22.2, -42.3 

7 18 Primary motor cortex Somato-motor R Parietal 40 52.8, —27.2, 40.9 
8 18 Premotor cortex Somato-motor R Parietal 1 32.4, —39.2, 49.6 
9 17 Anterior insula Cingular-opercular L Insula 13 —38.7, 8.1, -4.8 
10 16 Anterior insula Salience L Insula 45 —32.5, 22.1, 5.8 
11 16 Inferior frontal gyrus n/a L Prefrontal 11 —18.2, 12, -21 
12 16 Cerebellum n/a R Cerebellum n/a 6.1, —50.7, -12.3 
13 15 Cerebellum n/a L Cerebellum n/a —6.5, —50.2, -11.4 
14 14 Thalamus Subcortical L Subcortical n/a —4.9, —10.3, 5.8 
15 14 Thalamus Subcortical L Subcortical n/a —9.6, —25.4, —1.4 
16 14 Cerebellum n/a L Cerebellum n/a —10.3, —37.7, -25.1 
17 14 Hippocampus Default L Limbic n/a —32.1, —40.2, -4 
18 14 Retrosplenial cortex Default L Limbic 23 —8.6, —58.8, 17.6 
19 14 Thalamus Subcortical R Subcortical n/a —10.5, -26.8, -2.2 
20 14 Inferior frontal gyrus n/a R Prefrontal 11 15.6, 34.1, -22.6 
21 13 Brainstem Subcortical R Brainstem n/a —5, —21.5, -15.8 
22 13 Cerebellum n/a L Cerebellum n/a —21.2, —53.4, -23.6 
23 13 Midcingulate Default L Limbic 31 —8.8, —42.6, 50.1 
24 13 Fusiform gyrus Visual L Occipital 19 —25.9, —63.1, -12.3 
25 13 Brainstem n/a R Brainstem n/a 6.3, —24.9, -17.5 


BA, Brodmann area; K, degree; L, left; n/a, not available; R, right. 
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